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Abstract 

In the present work ferromagnetic ordering in the Hubbard model generalized by 
taking into account the inter-atomic exchange interaction and correlated hopping 
in partially filled narrow band is considered. In the case of weak electron-electron 
interaction the ground state energy and condition of ferromagnetic state realization 
are found by using the Green function technique. The obtained results indicate the 
important role of correlated hopping. 

PACS number(s): 71.10.Fd, 71.30.+h, 71.27.+a 

1 Introduction 

The problem of the origin of the metallic ferromagnetism, in spite of the variety of theo- 
retical attempts to solve it, still remains open. The simplest model at the basis of which 
one has tried to explain the physics of ferromagnetism in 3d transition metals is the 
Hubbard model that describes itinerant electrons in a single non-degenerate band inter- 
acting via on-site Coulomb repulsion |1J-||. The interplay between the kinetic energy 
of electrons and Coulomb interaction strongly influences the magnetic ordering. In some 
approximations the ferromagnetic solution in this model was obtained for large values of 
intra-atomic Coulomb interaction U ; in this connection note the Nagaoka's theorem [|| 
(the ferromagnetism in half-filled band with single electron or single hole at U — > 00). 
The problem of metallic ferromagnetism is studied intensively during last years (for recent 
reviews see Ref. || [|) however the undivided opinion is not reached so far. 

Nowadays we can distinguish a few ways to obtain the ferromagnetic solution: first 
is to stay in the framework of the Hubbard model but to use a special lattice geometry 
and density of states; second is to go beyond the Hubbard model including the band 
degeneracy (note that real ferromagnetic materials have the orbital degeneracy); and 
third is to take into account other matrix elements of electron-electron interactions in 
addition to the intra-atomic Coulomb interaction. We believe that third way can give the 
qualitative explanation of the magnetic properties of narrow band materials although the 
construction of consistent theory of metallic ferromagnetism should include the orbital 
degeneracy and real density of states. The aim of this article is to study the influence of 
so-called "off-diagonal" matrix elements of electron-electron interaction on the condition 
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of ferromagnetism stabilization. In spite of the fact that the magnitudes of the neglected 
in the Hubbard Hamiltonian terms are small in comparison with on-site repulsion or band 
hopping integral but they can essentially influence (assist or suppress) the ferromagnetic 
ordering. The importance of some of these matrix elements was pointed out in many works 
[[{]]- |f§; here we note the special role of the direct exchange interaction and correlated 
hopping (taking into account of the inter-atomic density-density Coulomb interaction 
which play essential role in charge ordering goes beyond the goal of this article). 

Last ten years the problem of importance of exchange interaction again is topical espe- 
cially in the works |I0|- ||13|| , where it was concluded that the nearest-neighbour exchange 
interaction J plays a fundamental role for the existence of ferromagnetism. In the cited 
papers the condition for ferromagnetism in the framework of mean-field theory has been 
obtained; also the exact diagonalization solutions for one-dimensional system with differ- 
ent band-fillings have been constructed and compared with results of mean-field theory 
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The exact criteria for the arising of ferromagnetism in the partial case of the model 



with correlated hopping have been derived by Strack and Vollhardt fL4| . These criteria 
also show the significance of the direct exchange interaction. Vollhardt and co-workers 
[p"5| , 16| within dynamical mean-field theory with the asymmetrical density of states having 



a peak near the band edge pointed out that already small direct exchange can stabilize 
the ferromagnetic ordering. 

As mentioned there is an additional mechanism which can favour ferromagnetic order- 
ing. The "off-diagonal" matrix elements of Coulomb interaction of electrons include the 
term that describes the density-dependent hopping - so-called correlated hopping. The 
importance of correlated hopping and its considerable place in the descriptions of metal- 
insulator transition and superconductivity are known ||18||- ||25|| . So we expect that taking 
into account the correlated hopping can give more clear understanding of the metallic 
ferromagnetism. The attempt to take into account the influence of correlated hopping 



has been done by Hirsch in Ref. ||13|| . It has been shown that correlated hopping can 
give rise to the asymmetry of the condition of ferromagnetic ordering (to enhance the 
tendency to ferromagnetism of electron-like versus hole-like carriers for small U, and have 
the opposite effect for large U). In our previous work we also have pointed out the impor- 
tance of electron-hole asymmetry in the theory of metal-insulator transition and metallic 
ferromagnetism [T^|. It is important that correlated hopping give rise to concentration 
dependency of hopping integral. 

The purpose of this paper is to investigate the possibility of ferromagnetic ordering 
and to study both the influence of correlated hopping of electrons and direct exchange 
interaction on the condition of ferromagnetism. The paper is organized as follows. In 
the Section 2 the Hamiltonian of the model with electron-hole asymmetry is formulated. 
In the Section 3 the single-particle Green function and energy spectrum are found in the 
case of weak electron-electron interaction. In the Section 4 the ground state energy of 
the model and the criteria of ferromagnetism stabilization are calculated and discussed. 
Section 5 is devoted to the conclusions. 
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2 The model Hamiltonian 



In accordance with the works |7|, [T7], [T?| we take into account in the Hamiltonian the 
matrix elements J(ijji) = J, J(ikjk) (k ^ i, k ^ j), J(iiij), where 
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J(ijkl) = j j 0*(r - R 4 )0(r - R fe )^-i— ^*( r ' - R,)0(r' - Rfidrdr . (1) 
The model Hamiltonian takes the following form |T^] 

H = - fi a ta a i° + Yl'ttv + nT i(^')) a i!r a ja + (T 2 (ij )a£ a j(T n id + /i.e.) 

ia ija ija 

+ U J2 n^nn + ' a ta a t*> a i°' a j^ ( 2 ) 

6 IJCFU 

where af a , a ia are creation and destruction operators of electron on site i, a =|, j, rii a = 
af a a ia , n = (n^ + n^), \i is the chemical potential, ty is the band hopping integral of 
an electron from site j to site i, U is the intra-atomic Coulomb repulsion and J is the 
exchange integral for the nearest neighbours, 

T 1 (i J ) = j2-mjk) 

k^i 

and T 2 (ij) = J(iiij) are the parameters of correlated hopping of electrons. The prime on 
the sums in Eq. (|2|) signifies that i ^ j. 

The peculiarities of the model described by the Hamiltonian (|2|) is taking into consid- 
eration the influence of the site occupation on the electron hoppings (correlated hopping), 
and the exchange integral. In this model an electron hopping from one site to another 
is correlated (in contrast to similar generalized Hubbard models) both by the occupation 
of the sites involved in the hopping process (with the hopping integral T 2 ) and the occu- 
pation of the nearest-neighbour sites (with the hopping integral T\) which we took into 
account by means of the Hartree-Fock approximation. 

Thereby the correlated hopping, firstly, renormalizes the initial hopping integral (it 
becomes concentration- and spin-dependent) and, secondly, leads to an independent on 
quasiimpulse shift of the band center, dependent on magnetic ordering. Taking into 
account the quantity of second order J is in principle necessary to describe ferromagnetism 
in this model. 

To characterize the value of the correlated hopping we introduce dimensionless param- 
eters Ti - Ml To - Ml 

exers t x - |ty| , r 2 - |ty| . 

3 The Green function and energy spectrum in case 
of weak interaction 

The single-particle Green function satisfies the equation 

(E - ii){{a pa \a+, a )) E = + J2 t ip( n )(M a i'cr))E 



3 



i i 

+2^T 2 (ip)((a^a i9 a pff \ap <T ))E + U((n p ^a pa \ap a )} E 

i 

+ J((a%> a p° ,a i°\ a t'°))E- ( 3 ) 

ja' 

Let us consider the system at weak intra-atomic Coulomb interaction (U < 2w, w = 
z\tij\ is the half of bare bandwidth, z is the number of nearest neighbours to a site). 
Then we can take into account the electron-electron interactions in the Hartree-Fock 
approximation: 

((a±ai«a ja \a+, a )) E ~ (a±a i& ) ((a ja \a+, a )) E ; (4) 
((a&ajjraurlap^E ^ (ata J - ff )<(a i(T |aJ ff ))B. 

We assume that averages {a~l a a ia ) = n a are independent of the number of the site (the 
uniform distribution of electronic density and magnetic moments is assumed). 

The approximation (|j) corresponds to the condition of weak electron-electron interac- 
tion, where doubly occupied states can be taken into account by the effective mean field 
depending on the electron concentration and correlated hopping integral. 

After the transition to Fourier representation we obtain for the Green function 

<M<v»k = ^ E _ l E ^y (5) 

where the energy spectrum is 

E a (k) = -fi + p a + n„U - zn a J + t(na)~f(k); (6) 

here the spin-dependent shift of the band center is 
2 

P° = ]y T 2 (y ) ( a t* a j°) > ( 7 ) 

ij 

7(k) = J2 e ,kR (the sum goes over the nearest neighbours to a site) and the spin and 
R 

concentration dependent hopping integral is 

t{na) = t + nT x + 2n- a T 2 + J^(a+,a jV >. (8) 

a' 

From the energy spectrum (0) by neglecting the last term in (Kp we obtain the result 
of papers (17|, [18| ; neglecting the correlated hopping Tiiij) - the result of Ref. |13 . 



The dependences of effective hopping integral on electron concentration and magneti- 
zation, a presence of the spin-dependent shift of band center are the essential distinctions 
of single-particle energy spectrum of the model described by Hamiltonian @) from the 
spectrum of the Hubbard model in the case of weak interaction. 

The concentration of electrons with spin o is 

+oo 

n a = ( p(e)f(e)de. (9) 
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Here p(e) is the density of states, /(e) is the Fermi distribution function. Let us assume 
the rectangular density of states: 

P( e ) = Jj E S(e - e(k)) = ^6(e 2 - w 2 ). (10) 

k 

In the case of zero temperature we obtain: 

e a + w 

here the value e a is the solution of the equation E c {e) = from which we obtain e a = 
where /x CT = fi — (3 a + zn a J — n^U and a a — 1 — Tin — 2r 2 n^ — ^ £}?v(l ~ n a>)- 

a' 

The system parameters are related by the equation 

t rr ft~ft n A 2 Jb(2 - n) - m 2 ] \ 

zJ + [/ + ^ — ^ = 2w 1 - Tin - r 2 1 v - ; *- . 12 

my 2w J 

The shift of band center is obtained from 
2 

ft = jyE^feK^) = -T2wna{n s - 1). (13) 
One can see that 

/3 i -/3 T = 2r 2 mu;(l-n). (14) 

In particular, the concentration dependent shift fll4T) of the centers of spin-up and spin- 
down electron bands at n < 1 is positive, at n > 1 is negative. In this manner we obtain 
the condition for the equilibrium value of magnetization: 

U z J 

1 - Tin - r 2 (2 -n) = — [1 + n(2 - n) - m 2 }. (15) 

2w 2w 

This result will be compared with result of next section. 

4 Ground state energy 

To calculate the ground state energy of the model we use the formula 

+oo 

£o = ^E / (*k(") + E)r{E)dE. (16) 

kcr -oo 

Here 

1 + exp 

is the spectral intensity of Green function (0). From Eq. (|HJ) one can obtain for the 
ground state energy the expression 

E o = 9 E {-^ria - (1 - Tin + a CT )n CT (l - 7v)w) , (18) 
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which can be rewritten in the form 

£o = 4 0) + 4 2) +4 4 \ (19) 
4 0) = | + §» + - (1 - + r 2 )n)(2 - , 

4 2) = (l - rm - r 2 (2 - n) - ^ - ^(1 + n(2 - n))) f m 2 , 
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One can see that the condition for equilibrium value of magnetization obtained from 
expression dE Q /dm = is equivalent to Eq. (plj). The position of the minimum of ground 
state energy depends on values of model parameters. In Fig.l the energy difference be- 
tween paramagnetic and ferromagnetic states is plotted as a function of the magnetization. 
From this plot we can see that at some values of the parameters a ferromagnetic ordering 
occurs with increase of inter-atomic exchange integral. The values of magnetization at 
which the ground state energy of the system has a minimum depends on magnitude of 
inter-atomic exchange. The condition of ferromagnetic ordering stability d 2 Eo/dm 2 < 
can be obtained as 

zJ_ l-nr 1 -r 2 (2-n)-£ 

2w l + n(2-n) K } 

For the case of zJ/2w = from Eq. (|20|) we obtain a generalization of the Stoner criterion 
which takes into account the correlated hopping 

Up{e F ) > 1 -Tin -T 2 (2-n). (21) 

From the condition of the minimum of ground state energy dEo/dm = one can also 
obtain for the value of exchange integral at fixed magnetization m 



zJ 



zJ 1 -nri - r 2 (2 -n) - 



',22) 



2w 1 + n(2 — n) — m 2 

The condition of full spin polarization (m = n) is 

zJ_ l-nr 1 -r 2 (2-n)-£ 
2w l + 2n(l-n) ' 

From Eq. (E0)-(P3) one can conclude: the competition between correlated hopping pa- 
rameters Ti and r 2 determines more favourable situation for the ferromagnetic ordering. 
For example, if r 2 > T\ then the systems with the electron concentration n < 1 are more 
favourable to ferromagnetism than the systems with n > 1 and vice versa (see Fig. 2, 3; 
the values above the line correspond to the ferromagnetic state of the system, below the 
line - to the paramagnetic one), but the correlated hopping parameter r 2 influences the 
condition of ferromagnetic ordering stronger than the parameter t\ (see Fig. 4). Note 
also that the region of partial spin polarization is narrowed with a deviation of electron 
concentration from half-filling (see Fig. 2, 5, 6). 

One can see that in the less than half-filled band, correlated hopping leads to the 
stabilization of ferromagnetism as well as inter-atomic exchange interaction and intra- 
atomic Coulomb interaction; the larger is the electron concentration n, the smaller is the 
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critical value of the exchange integral for the occurrence of ferromagnetism. These our 



results are in accordance with the results of Ref. [[L3|] . 

The magnetization found from the condition of minimum of the ground state energy 
is written as 

m = d)" + n(2 - "» + s; + nn + M2 - n) - ' • (24) 



The magnetization defined by Eq. (j24|) is plotted as a function of zJ/w for different values 
of U/w in the Fig. 7. In the ferromagnetic state the magnetization rapidly increase with 
zJ/w and reach its maximum value at relatively small values of exchange integral what 
indicates the importance of inter-atomic exchange interaction for ferromagnetism of the 
system. The dependence of the magnetization on concentration of electron is shown in 
Fig. 8. At fixed values of model parameters the change of electron concentration (doping 
of the system) can induce ferromagnetism. 

In Fig. 9 the dependence of paramagnetic and ferromagnetic ground state energies of 
the system on zJ/w is plotted for different values of n. For all values of the electron 
concentration at sufficiently large values of zJ/w the energy of the ferromagnetic state 
lies much lower then the paramagnetic one. The position of the transition point depends 
on values of U and r 1; r 2 . With the increase of these parameters the critical value of inter- 
atomic exchange interaction becomes smaller (the influence of the Coulomb repulsion is 
illustrated in Fig. 10). 

The energy difference between the paramagnetic and ferromagnetic states, the value 
of magnetization as functions of band filling are plotted in Fig. 11. Depending on the 
value of n the state of the system can be para- or ferromagnetic. For the zero values of 
correlated hopping parameters (the curve 1 on Fig. 11 (a) ) in the system with increase 
of n the ferromagnetic order occurs. Note that with further increase of n the system 
becomes fully spin-polarized (see fig. 11(b) ) however in almost half-filled band the spin 
polarization is partial again. In the more then half-filled band (in consequence of the 
electron-hole symmetry of this case) the noted magnetic behaviour repeats in inverse 
sequence. In the case of non-zero values of correlated hopping parameters (the curves 
2,3,4 on Fig. 11 (a) ) we have the valuable decrease of ferromagnetic ground state energy 
comparing to paramagnetic one; the influence of correlated hopping t\ is stronger at n > 1 
and the influence of r 2 is stronger at n < 1. In the case of equal values of T\ and r 2 the 
electron- hole symmetry retrieves (the curve 4 on Fig. 11 (a)). 



5 Conclusions 

Although the consistent theory of ferromagnetic ordering in transition metal compounds 
can be constructed only in the model including the orbital degeneracy of the band, the 
qualitative character of observed properties can be interpreted in the framework of gen- 
eralized model of non-degenerate band. We have studied the ferromagnetism in this 
model of at weak intra-atomic interaction. The peculiarity of the studied model is taking 
into account besides the intra-atomic Coulomb and inter-atomic exchange interactions 
the additional terms describing the interaction of the hopping electrons with other elec- 
trons (correlated hopping). In the considered case of weak electron-electron interaction 
the inter-atomic exchange interaction leads to the transition from a paramagnetic to a 
ferromagnetic state of the system. 
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The critical value of exchange interaction strongly depends on model parameters. The 
correlated hopping causes the spin-dependent shift of the band center which can lead to 



the ferromagnetism. This our result is in accordance with the results of Refs. [|iT| , [H| 



In distinction from the model considered in Ref. Jl3| our Hamiltonian contains additional 
mechanism of correlated hopping which is shown to favour the ferromagnetic ordering. 
The intra-atomic Coulomb interaction is the factor favouring the ferromagnetic ordering 
but this interaction itself can not stabilize the ferromagnetic ordering if U < 2w (see 
Fig.l). The ground state energy has strong filling dependence and its behaviour can be 
changed by rise of the correlated hopping parameters. 

To conclude, in the case of weak electron-electron interaction the inter- atomic exchange 
interaction plays a substantial role in stabilization of the ferromagnetic ordering; the intra- 
atomic Coulomb interaction and the correlated hopping also favour ferromagnetism. The 
non-equivalence of less then half-filled and more than half-filled cases (concerning the 
condition of ferromagnetism realization) is shown which is characteristic for the models 
with correlated hopping of electrons. 
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Figure captions 

Fig.l The energy difference between paramagnetic and ferromagnetic states as a 
function of the magnetization m at ^ = 1, T\ = r 2 = 0.2 and n = 0.7. Upper, middle and 
lower curves correspond to — = 0, — = 0.3 and — = 0.5 respectively. 

1 w ' w w 

Fig. 2 Ground state phase boundaries at — = 0.8, T\ = 0.15 and r 2 = 0.2: curve 1 
corresponds to the paramagnetic case, curves 2,3 - the cases of partially spin polarization, 
50 and 75 percents respectively, curve 4 corresponds to the case of full spin polarization. 

Fig. 3 The critical value of the exchange integral as a function of r = r 2 at ^ = 0.8 
and T\ = 0.15: curve 1 corresponds to n — 0.5, curve 2 corresponds to n — 1 and curve 3 
corresponds to n — 1.5. 

Fig. 4 The critical value of the exchange integral as a function of electron concentration 
n for given magnetization at T\ = 0.15, r 2 = 0.2 (solid lines) and T\ = 0.2, r 2 = 0.15 
(dashed lines), ^ = 0.8: curves 1,2,3 and 4 correspond to m — 0.1, m = 0.5, m = 0.7 and 
m = 0.9 respectively. 

Fig. 5 Lines of constant magnetization for half-filled band: curves 1,2 and 3 correspond 
to the cases m = 0, m = 0.7 and m — 1 respectively. 

Fig. 6 Lines of constant magnetization for quarter-filled band: curves 1,2 and 3 corre- 
spond to the cases m = 0, m = 0.7 and m = 1 respectively. 

Fig. 7 The magnetization m as a function of — at T\ = r 2 = 0.1 and n = 0.8. Curves 
1, 2 correspond to — = 1.2 and — = 0.7 respectively. 

Fig. 8 The magnetization m as a function of n at — = 1.2 and t± — r 2 = 0.1. Solid 
curve corresponds to — = 0.3, dashed curve corresponds to — = 0.6. 

Fig. 9 The paramagnetic (dashed line) and ferromagnetic (solid line) ground state 
energies of the system as a function of ^ at ^ = 0.5, n = t 2 = 0.3. Upper curves 
correspond to n — 1.2, lower ones correspond to n — 0.7. 

Fig. 10 The ground state energy as a function of — at T\ = r 2 = 0.2 and n = 0.8. 
Curves 1,2 and 3 correspond to — = 1, — = 0.6 and — = 0.2 respectively. 

Fig. 11 The energy difference between the paramagnetic and ferromagnetic states 
(a) and the value of magnetization m (b) as functions of band filling n at ^ = 1.2, 
^ = 0.45. Curves 1, 2, 3 and 4 correspond to the cases T\ = r 2 = 0, T\ = 0, r 2 = 0.05, 
Ti = 0.05, r 2 = and T\ = r 2 = 0.05 respectively. 
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